Collinear helium under periodic driving: stabilization of the asymmetric stretch orbit 
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The collinear eZe configuration of helium, with the electrons on opposite sides of the nucleus, is 
studied in the presence of an external electromagnetic (laser or microwave) field. We show that the 
classically unstable "asymmetric stretch" orbit, on which doubly excited intrashell states of helium 
with maximum interelectronic angle are anchored, can be stabilized by means of a resonant driving 
where the frequency of the electromagnetic field equals the frequency of Kepler-like oscillations along 
the orbit. A static magnetic field, oriented parallel to the oscillating electric field of the driving, 
can be used to enforce the stability of the configuration with respect to deviations from collinearity. 
Quantum Floquet calculations within a collinear model of the driven two-electron atom reveal the 
existence of nondispersive wave packets localized on the stabilized asymmetric stretch orbit, for 
double excitations corresponding to principal quantum numbers of the order of N-Z 10. 





I. INTRODUCTION 

The correlated dynamics of two-electron atoms under 
external electromagnetic driving represents a fascinating 
topic of atomic physics. The outstanding example in this 
context is double ionization of helium in the presence of a 
strong laser pulse. At not too high intensities of the laser 
(near 10 15 W/cm 2 [10), the electrons are, contrary to 
natural expectation, not emitted in a sequential process 
which was clearly manifest in the original experiments 
111 Q- Further experimental investigations (e.g 
and a number of theoretical studies (e.g., 
I) have revealed that a nontrivial multistep process 
is responsible for nonsequential double ionization in this 
intensity regime, involving tunnel ionization of one of the 
electrons, recollision of this outer electron with the ionic 
core, followed by subsequent excitation and ionization of 
the inner electron. 

In the regime of highly excited states, complex dy- 
namics can also take place in the presence of moderate 
electromagnetic fields, which are not strong enough to 
fully ionize the atom, but substantially perturb the elec- 
tronic motion. A striking example in this context is the 
occurence of nondispersive wave packets in microwave- 
driven hydrogen: Rydberg wave packets, which are nor- 
mally subject to spreading and collapse after a limited 
number of Kepler cycles, can be stabilized and "kept in 
shape" over a practically arbitrary amount of time by 
means of a resonant electrical driving (with linear pdlll2j 
or circular polarization 0,0,^1). This field needs to 
be phase-matched with the wave packet in such a way 
that components associated with deviating classical tra- 
jectories are driven back to the resonant Kepler orbit. 
For a wave packet that performs parabolic-state-like os- 
cillations along the axis of a linearly polarized field for 
instance 0, 0] > this means that the external force on 
the electron needs to be directed towards the nucleus at 
the outer turning point of the classical motion, and away 
from it at the inner turning point (i.e., when the elec- 



tron collides with the nucleus). The wave packet then 
corresponds to a Floquet eigenstate of the periodically 
driven system, which is localized in phase space on the 
regular island that is associated with the above nonlin- 
ear resonance. The coupling to the "chaotic sea" outside 
the island, and subsequently also to the ionization con- 
tinuum, is mediated via a classically forbidden tunneling 
process through the dynamical phase space barriers of 
the island, what leads to almost "eternal" lifetimes of 
such wave packets [H[H0. 

The concept of nondispersive wave packets, which is of 
potential interest in the context of quantum control, can 
be generalized to the correlated dynamics in two-electron 
atoms. This was indeed shown for the classically stable 
"frozen planet" configuration of helium 0, 0], where 
both electrons are located on the same side of the nu- 
cleus with asymmetric excitation. The application of a 
linearly polarized microwave perturbation to this config- 
uration induces analogous islands of regular motion in 
phase space, which arise here from a resonant driving of 
the outer electron 0,0]. A static electric field, how- 
ever, is needed to ensure classical stability with respect 
to deviations from collinearity. The existence of nondis- 
persive two-electron wave packets localized on these res- 
onance islands has indeed been revealed in Floquet cal- 
culations within a quantum model that represents the 
analog of the collinear two-electron configuration [20ll2l| . 

In this contribution, we examine to which extent the 
mechanism that leads to nondispersive wave packets can 
be used to stabilize a doubly excited configuration - i.e., 
to enhance the lifetime of the associated autoionizing 
states by increasing the amplitude of the external elec- 
tromagnetic field. Originally, the phenomenon of sta- 
bilization was discussed in the context of one-electron 
atoms in ultra-intense high-frequency laser fields ("adi- 
abatic stabilization" ) , where it was shown that the ion- 
ization rate of the atom decreases with increasing laser 
intensity |23 (without, however, taking into account rel- 
ativistic effects; see 23]). The effect that we are aiming 
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at is more of intermediate nature (i.e., it manifests itself 
at moderate intensities of the electromagnetic field) and 
relies on the structure of the underlying classical phase 
space. The basic question is to which extent the driving 
field can create dynamical phase space barriers (invari- 
ant tori) around a periodic orbit that would be unstable 
without the external perturbation. Floquet eigenstates 
that are locally quantized on this stabilized orbit are then 
semiclassically protected against decay, which enhances 
their lifetime compared to the unperturbed atom. 

A natural candidate for this stabilization mechanism 
is the "asymmetric stretch" orbit [24|. In this orbit, the 
electrons are located on opposite sides of the nucleus and 
perform slightly perturbed Kepler oscillations with oppo- 
site phase; whenever one of the electrons hits the nucleus, 
the other one reaches the outer turning point of the orbit. 
The dynamics along this orbit is stable with respect to 
transverse bending perturbations, but unstable against 
deviations in axial direction. As was shown by Richter 
and Wintgen |25j . intrashell resonances with maximum 
interelectronic angle are predominantly scarred along this 
orbit - which can be seen as a consequence of the promi- 
nent role that this orbit plays in the semiclassical quan- 
tization of helium [24[ • 

Due to the phase shift between the Kepler-like oscilla- 
tions, the stabilization mechanism that was used to cre- 
ate nondispersive wave packets in driven hydrogen can 
now be applied simultaneously to both electrons: We ex- 
pose the collinear configuration to a coaxially polarized 
electromagnetic field the force of which is directed to- 
wards the nucleus at the outer turning point of each 
electron, and outwards when the electron collides with 
the nucleus. Indeed, we shall see in Section [H] that the 
asymmetric stretch orbit can thereby be stabilized within 
the phase space of collinear motion, if the amplitude of 
the driving field is appropriately chosen. In the presence 
of the time-periodic driving, the orbit is no longer sta- 
ble with respect to deviations from collinearity - due to 
mixing of states with different angular momentum. We 
shall argue in Section II I II however, that a static mag- 
netic field can be used to enforce the transverse stabil- 
ity of the driven configuration. In Section I1VI finally, 
the basic properties of nondispersive two-electron wave 
packets localized on the stabilized asymmetric stretch or- 
bit are elaborated using a one-dimensional model that 
represents the quantum analog of the collinear eZe con- 
figuration. We shall point out that these wave packet 
states evolve diabatically from the unperturbed asym- 
metric stretch state, and that their lifetimes exhibit a 
local maximum at finite values of the field amplitude. 



II. CLASSICAL STABILIZATION OF THE 
COLLINEAR ASYMMETRIC STRETCH ORBIT 

In atomic units, which are used throughout this pa- 
per, the classical Hamiltonian of the electromagnetically 



driven helium atom reads 

" 2 + 2 |n| |r 2 | + In-ral 

+ Oi + zi)Fcos{ut + ip). (1) 

Here, rj = {x l ,y l ,z i ) and p t = (p X i,P y i,Pzi) denote the 
position and momentum of electron i — 1,2, respectively, 
Z = 2 is the nuclear charge, and F and u> represent 
the amplitude and frequency of the external driving field 
which is linearly polarized along the z axis. In analogy 
to driven hydrogen 26] as well as to the unperturbed 
helium atom |27j . the Hamiltonian exhibits general 
scaling laws: The classical dynamics generated by (QJ 
remains invariant if all variables and parameters of the 
system are transformed according to 



r, .— > v 2 v l (z = 1,2), (2a) 

Pl .— > I/" 1 Pi (z = l,2), (2b) 

t i — > v 3 t, (2c) 

F i — ► v~ 4 F, (2d) 

uj i — ► v^uj, (2e) 

H i — ► v~ 2 H, (2f) 



where v represents an arbitrary, real positive quantity. 
We shall therefore restrict the classical analysis to a fixed 
value w = 1 of the driving frequency, and use the above 
scale invariance J2J to deduce the actual dynamics at the 
energy range of interest. Fq denotes, in the following, 
the value of the "scaled" field amplitude F at driving 
frequency lu = 1. 

In the following, we shall concentrate on the invariant 
subspace of collinear motion along the z axis, with the 
electrons on opposite sides of the nucleus. Such a collinear 
eZe motion is generated by initial conditions of the form 
Xi(t = 0) = y t {t = 0) = Pxt (t = 0) = p yi (t = 0) = for 
i = 1,2, and z%(t = 0) > whereas z%(t = 0) < 0. Due to 
the Coulomb singularity in the interaction, the electrons 
cannot pass the nucleus at the origin. This implies that 
zi(t) > and z 2 (t) < for all times t > 0. 

To ensure a stable numerical integration of the classical 
equations of motion, we perform a Kustaanheimo-Stiefel 
transformation 28]. This procedure, which is described 
in Appendix ^ for the special case of collinear motion, 
introduces new phase space variables which do not di- 
verge at electron-nucleus collisions. Triple collisions - 
i.e., the simultaneous encounter of both electrons at the 
nucleus - cannot be regularized in such a way; they repre- 
sent "true" singularities of the dynamics, where different 
manifolds of trajectories merge together without a well- 
defined continuation [2j|. 

Correspondingly, the Kolmogorov-Arnold-Moser 
(KAM) theorem cannot be applied to this system, 
which is reflected by the fact that the classical dynamics 
of the unperturbed three-body problem (i.e., in the 
field-free case Fq = 0) is fully chaotic, even in the limit 
of very weak electron-electron interaction (which would 
correspond here to very small l/Z) 0,113. Indeed, it 
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FIG. 1: a) The asymmetric stretch orbit of the collinear eZe 
configuration. Plotted are the coordinates of the two elec- 
trons (zi > and 22 < 0) as a function of time, with the 
nucleus resting at z = 0. The scaling (see Eq. (|5J) is chosen 
such that the frequency of the orbit equals u = 1. 
b) A small change of the initial condition (the "upper" elec- 
tron is displaced by Azi = 0.01) leads to an appreciable devi- 
ation from the asymmetric stretch orbit (shown in grey in 
the background) within a couple of oscillation cycles. At 
t/(2n) ~ 6, the configuration ionizes after a nearby triple 
collision. 



is possible to represent each periodic orbit by a unique 
sequence of symbols, which is basically determined by 
the order in which the electrons hit the nucleus |24| . 
The simplest orbit with respect to this symbolic code 
is the "asymmetric stretch" orbit, shown in Fig. ^a), 
in which the electrons collide with the nucleus in a 
perfectly alternating way. Since in this configuration the 
electrons avoid a simultaneous encounter at the nucleus 
in the best possible way, the stability exponent of this 
orbit is rather low compared to other, more complicated 
periodic orbits |24j . 

Nevertheless, a small deviation from the orbit's ini- 
tial condition leads to disintegration of the configura- 
tion on rather short time scales. This is illustrated 
in Fig. ^b), which shows what happens when one of 
the electrons is displaced in configuration space by an 
amount of Az = 0.01: After a few cycles, the Kepler-like 
oscillations of the electrons get out of phase, which leads 
to a nearby triple collision where a large amount of energy 
is transferred between the electrons, and to subsequent 
(single) ionization of the atom. 

The external time-dependent electric field is now ap- 
plied to this configuration in order to compensate this 
destabilization phenomenon. As in the case of nondisper- 
sive wave packets in driven hydrogen, the relative phase 
between the field and the electronic motion is adjusted 
such that the field forces the electron towards the nucleus 
at the outer turning point of the Kepler-like oscillation, 
and away from it at the inner turning point. In this way, 
deviations from the periodic orbit are counterbalanced 




FIG. 2: The "Lyapunov exponent" a of the asymmetric 
stretch orbit, i.e. the logarithm of the largest eigenvalue of 
the stability matrix, is plotted as a function of the field am- 
plitude F . In between Fo ~ 0.044 and F ~ 0.154, a is 
negative, which indicates that the asymmetric stretch orbit 
is stable and constitutes the center of a regular island in the 
phase space of driven helium. 



by the field. Since their oscillations take place with op- 
posite phase, the phase-matching condition can be simul- 
taneously satisfied for both electrons. For the particular 
type of initial condition in Fig. ^ where zi(0) = —22(0) 
and Pzi(0) = Pz2(0) < 0, we need to choose <p — ir/2 as 
initial field phase in Eq. Q). 

Due to the high dimensionality of the extended phase 
space, which is spanned by the positions and momenta 
of the electrons and by the phase of the external driving, 
regular islands that result from this stabilization mecha- 
nism cannot be visualized by means of Poincare surfaces 
of section. To identify the periodic orbit at finite field am- 
plitude Fo 7^ 0, we employ an iterative procedure which 
is described in Appendix [5] The key ingredient to this 
procedure is the fact that initial conditions of asymmetric 
stretch orbits at different field amplitudes Fq and Fq+5Fq 
lie close to each other in phase space - i.e., within each 
other's linear neighbourhood - as long as \SFq\ is rather 
small. Hence, by varying Fo in sufficiently small steps 
and using a Newton- Raphson method to adapt the initial 
condition, we can "trace" the periodic orbit as a function 
of the field amplitude. 

Fig.|21shows the Lyapunov exponent a = In X ma x of the 
driven asymmetric stretch orbit, which is calculated from 
the largest eigenvalue \ ma x of the stability matrix associ- 
ated with one oscillation period. In between Fq ~ 0.044 
and Fq ~ 0.154, this eigenvalue is smaller than unity, 
which formally results in a negative Lyapunov exponent 
and indicates stable dynamics in the vicinity of the orbit. 
This is indeed confirmed by plotting the corresponding 
trajectories in configuration space - e.g., for Fo = 0.074 
as shown in Fig. |3| Instead of increasing exponentially 
with time, a small deviation from the periodic orbit leads 
to stable, quasiperiodic oscillations around the orbit. 

A semiclassical estimation of the minimum atomic 
excitation at which this regular island supports fully 
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FIG. 3: a) The stabilized asymmetric stretch orbit at Fo = 
0.074 and u = 1, plotted over 20 field cycles. Initial condition: 
Zi = -z% = 2.0367, p z i = p 22 = -0.54512. Plotted in grey in 
the background is a trajectory with a slightly different initial 
condition (Azi = 0.1). The effect of such a deviation from the 
fundamental periodic orbit is magnified in (b) in the vicinity 
of the upper electron's outer turning point (on top of the 
fundamental orbit which is plotted in grey): In contrast to 
the field-free case (Fig.0, the electrons maintain their relative 
phase and perform stable oscillations around the asymmetric 
stretch orbit. 



localized quantum states is provided by the Einstein- 
Brillouin-Keller (EBK) quantization criterion (see, e.g., 
IHEl). The latter states that any quantized torus must 
fulfill 



ptfq = 27iF ( r 



(3) 



for all closed curves C that are contained within the sur- 
face of the torus. Here, (p, q) are canonically conjugate 
phase space variables, n > is an integer, and p repre- 
sents the Maslov index accumulated along C. If we choose 
the curve C to encircle the interior of the torus within the 
subspace that is spanned by a particular pair (p.; , qi ) of 
canonically conjugate variables (which implies \i = 2), 
we obtain the requirement that the "cross section" area 
Aj = § Pidqi enclosed by the torus within that subspace 
must equal 2irh(n + 1/2) for some n > 0. Hence, in order 
to obtain at least one quantized state within the island, 
we need to require that its cross section area with respect 
to any pair (p.; , qA of phase space variables is at least of 
the order of ttH; in that case, at least one of the invariant 
tori would fulfill the EBK criterion for n = 0. 

Fig. 01 visualizes the cross section through the asym- 
metric stretch island at Fo = 0.074. From a discrete 
lattice of 2867 equidistant initial conditions within the 
range 1.32 < z\ < 2.52 and 0.24 < p zl < 1.16 (the initial 
values of z^, and p Z 2 are fixed to the ones associated with 
the fundamental periodic orbit), only those are plotted 
that lead to stable quasiperiodic motion after 1000 field 
cycles. From the number n s = 447 of such stable ini- 
tial points, we obtain the cross section area of the island 
within the Z\—p z \ space (and, by symmetry, also within 



FIG. 4: Cross section through the regular island around the 
stabilized asymmetric stretch orbit at Fo = 0.074. Plotted 
are initial conditions in the Z\ - p z \ subspace that lead to 
stable quasiperiodic motion after 1000 field cycles, while the 
initial condition of the other electron is fixed to the one of the 
fundamental periodic orbit: 22 = 2.0367, p Z 2 = —0.54512. z\ 
and Pzi are varied in discrete steps of the size Az\ = Ap z i = 
0.02. The cross section area of the island extracted from this 
figure equals Aq — 0.179. 
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FIG. 5: Cross section area Aq of the stabilized asymmetric 
stretch island as a function of the field amplitude Fo, calcu- 
lated in the same way as in Fig. [4] The global maximum is 
encountered at Fq = 0.074. 



the Z2~Pz2 space) according to 

A = n s AziAp zl ~ 0.179 



(4) 



with Azi = Ap z \ = 0.02 the spacing between adjacent 
lattice points in position and momentum, respectively. 
This is indeed the largest cross section area of the is- 
land that can be obtained by stabilizing the asymmetric 
stretch orbit. Fig.[5]shows the cross section area Aq, cal- 
culated in the above way, as a function of the field ampli- 
tude Fo; clearly, the maximum occurs around Fq = 0.074. 

The integrated action along the asymmetric stretch or- 
bit is given by 



S 



J a 



i=l,2 ' : ' 



(5) 



From numerical integration, we obtain S ~ 19.0 = Sq at 
uj = 1. If we assume (in analogy to the quantization of 
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Kepler orbits in hydrogen) that semiclassically S ~ AirN 
for intra-shell asymmetric stretch states with principal 
quantum number N (h is dropped again here and in the 
following) , and if we take into account that phase space 
cross section areas § pdq scale in the same way as action 
integrals when a scaling transformation according to (J2J) 
is performed, we obtain that a cross section area Aq at 
the scaling corresponding to u> = 1 is equivalent to a cross 
section area A — 4:TtNA / Sq at the scaling corresponding 
to the state with principal quantum number N. Hence, 
requiring that A > 7r, we formally obtain the minimum 
principal quantum number 

N min = S /(4Ao) ~ 27 (6) 

at which fully localized quantum states are to be expected 
on the stabilized asymmetric stretch island. 

It should be noted that the above estimation JJjJ is 
fairly imprecise and has to be taken as approximate 
guideline for the order of magnitude of N m i n , rather than 
as a precise criterion. This is, on the one hand, due to 
the simplifications that are involved in determining the 
cross section area of the outermost invariant torus of the 
island. On the other hand, the notion of an "outermost" 
torus itself is, strictly speaking, not meaningful in five or 
more dimensions, where invariant tori (with codimension 
larger than unity) do not divide the phase space in dis- 
connected segments. As a consequence, chaotic sublayers 
"within" the island are connected with each other and 
with the chaotic sea "outside" the island, and trajecto- 
ries starting on such sublayers can, via Arnold diffusion 
[32j| . leave the island on finite (though typically rather 
long) time scales. This is illustrated in Fig. EJ which 
shows the time evolution of a trajectory starting close 
to the boundary of the island. After a seemingly stable 
quasi-periodic oscillation over more than 80 field cycles, 
the correlation between the electrons breaks down and 
the atom ionizes. 

At finite values of h, moreover, the presence of Cantori 
|33j can inhibit quantum transport also in the immediate 
vicinity outside the regular island (where "hierarchical 
states" are localized |2J|)- The effective size of the local- 
izing region that the quantum system "sees" may there- 
fore be considerably larger than the size of the classical 
island. Correspondingly, estimations via the EBK crite- 
rion are typically rather conservative and tend to overes- 
timate the actual value of the minimum scaling needed 
to obtain a localized state on the island (an extreme case 
was reported in |35l l3r| where quantum states associated 
with the "Langmuir" orbit of helium were semiclassically 
predicted for N > 500 and quantum mechanically found 
at N = 10). Indeed, we shall see in Section HV1 that 
double excitations with principal quantum numbers of 
the order of N ~ 10, which are experimentally accessible 
[37i Isli ] , would be sufficient to obtain nondispersive wave 
packets anchored on the asymmetric stretch orbit. 
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FIG. 6: Long-time instability of a trajectory starting close to 
the boundary of the regular island at Fq = 0.074. The initial 
conditions are as in Fig. EI a), with the upper electron dis- 
placed by Azi = 0.2. After a seemingly stable quasi-periodic 
oscillation over about 85 field cycles, the correlation between 
the electrons "suddenly" breaks down. A nearby triple col- 
lision transfers then a large amount of energy between the 
electrons, which leads to ionization. 

III. DEVIATIONS FROM COLLINEARITY 

In contrast to deviations within the collinear subspace, 
the transverse degrees of freedom do not lead to destabi- 
lization of the unperturbed asymmetric stretch orbit. If 
both electron are symmetrically displaced in a direction 
that is perpendicular to the z axis - e.g., by nonvanishing 
initial components y\{t = 0) = y^it = 0)) - the effective 
repulsion between them is enhanced due to the reduced 
screening by the nucleus, which drives the electrons back 
to the z axis and leads to stable bending vibrations of the 
configuration (see in this context also the "asynchronous" 
orbit discussed in |3jJ, |40( ) . Marginal stability (with Lya- 
punov exponent 0) is encountered for antisymmetric dis- 
placements - e.g., with yi(t = 0) = —yi{t = 0) - which 
is a consequence of angular symmetry and the resulting 
conservation of the total angular momentum. 

In the presence of an external driving, however, with 
the phase of the field chosen such as to stabilize the 
collinear orbit, the angular symmetry is broken and the 
marginally stable degree of freedom is transformed into 
an unstable one. This can be understood from the net 
torque that the field exerts onto the atom. If the config- 
uration is slightly tilted with respect to the z axis (e.g., 
by small initial components yi(t = 0) = — y^{t = 0)), 
the external force is no longer directed parallel to the 
semimajor axis of the Kepler orbits, but exhibits a small 
perpendicular component. At the outer turning point of 
the Kepler oscillation, where the contribution to the net 
torque is most effective, this perpendicular component of 
the force is directed away from the z axis (see Fig. Ufa)) 
and therefore tends to enhance the initially small dis- 
placement from collinearity. 

This phenomenon of transverse instability arises also 
in the resonantly driven hydrogen atom. In the one- 
electron case, however, the two-dimensional motion that 
results from a small displacement from the field polariza- 
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FIG. 7: Transverse instability of the driven asymmetric 
stretch orbit. Plotted is, in the y-z plane, the time evolu- 
tion of a trajectory with initial condition in the vicinity of 
the collinear stabilized orbit at Fo = 0.074 (the same initial 
conditions as in Fig.[3Ja), with the electron on the right-hand 
side transversally displaced by Ay = 0.01). The trajectory 
is shown in the time intervals (a) - 5 field cycles, (b) - 
28 field cycles, (c) 28 - 38 field cycles, and (d) 38 - 50 field 
cycles. 



tion axis is bounded by means of dynamical barriers in 
phase space and does not lead to ionization. Such dynam- 
ical barriers do not exist in the two-electron case, due to 
the high dimensionality of the phase space. Fig. {7\ shows 
what happens when the configuration starts from an ini- 
tial condition that corresponds to the stabilized collinear 
orbit, with one of the electrons slightly displaced in y 
direction. We clearly see that the external field drives 
the configuration away from the z axis (Fig. 0a)) and 
rotates it around the nucleus (Fig. Efb,c)). At about 40 
field cycles, when one such rotation is completed, the cor- 
relation between the electrons breaks down and the atom 
eventually undergoes single ionization. 

The scenario encountered here is indeed very similar 
to the case of the resonantly driven Zee (frozen planet) 
configuration [iH IT^. l2l| , which is also unstable against 
deviations from collinearity. For that case as well as 
for the resonantly driven hydrogen atom |4lL |42| . it is 
known that the application of an additional, static elec- 
tric field parallel to the polarization of the driving can 
enforce transverse stability of the configuration. For this 
purpose, the static electric field needs to be directed such 
that it prevents the electrons from approaching the z = 
plane - and thereby counterbalances the destabilization 
mechanism induced by the torque of the time-periodic 
driving. 



0.6 



y 




FIG. 8: Transverse stabilization of the asymmetric stretch 
orbit at Fo = 0.074 by means of a magnetic field B — 0.5 
polarized along the z axis. The initial condition corresponds 
to the stabilized collinear orbit, with the electron on the right- 
hand side transversally displaced by Ay = 0.5. Plotted are 
the y and z components of the three-dimensional trajectory 
over 30 field cycles. We see that the transverse confinement 
provided by the magnetic field leads to a stable, quasiperiodic 
oscillation around the field polarization axis. 

In the case of the driven eZe configuration, a homo- 
geneous static electric field cannot be used to prevent 
transverse destabilizaton, since it cannot be oriented in 
such a way that it simultaneously keeps both electrons 
from approaching the z = plane. An effective alterna- 
tive, however, is provided by the possibiliy of adding a 
strong static magnetic field, which is oriented along the z 
axis, i.e., parallel to the polarization of the driving. Due 
to the tight diamagnetic confinement created by such a 
magnetic field, small deviations from collinearity are not 
enhanced, but lead to stable (three-dimensional) rota- 
tions of the electrons around the z axis. This is illustrated 
in Fig. |21 which shows the two-dimensional projection of 
a trajectory that starts in the vicinity of the collinear 
stabilized asymmetric stretch orbit at Fo = 0.074, in the 
presence of a magnetic field with strength Bq = 0.5 ori- 
ented along the z axis. The initial displacement of the 
electron on the right-hand side is y\ = 0.5, which indi- 
cates a rather large phase space volume of the associated 
regular island. 

For the asymmetric stretch orbit at Fo = 0.074, we find 
that magnetic field strengths in the range Bo ~ 0.2 . . . 2 
can be used to provide transverse stability of the config- 
uration |43[ . Since a homogeneous magnetic field B is 
incorporated into the Hamiltonian via the substitution 

Pi^Pi + -(Bxr) (7) 

(i = 1,2), it is scaled according to B 1 — > v~ 3 B 
when a scaling transformation of the type (J2J) is ap- 
plied. Hence, at double excitations corresponding to 
the principal quantum number N = 25 . . . 30, where the 
EBK criterion predicts fully localized quantum states on 
the asymmetric stretch island, a "scaled" magnetic field 
strength B = 0.2 (at uo — 1) would roughly correspond 
to B = (4ttN/S q )- 3 B ~ 10 Tesla (with S ~ 19, see 
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Section [Hj. Similarly, the resonant driving frequency 
would be given by ui = {4ttN/S )~ 3 ~ 2tt x 10 12 Hz. 
To realize the field amplitude F = (AttN/ S Q )- 4 F with 
Fo ~ 0.07, radiation intensities of the order of I ~ 10 5 
W/cm 2 would be required. This points towards quantum 
cascade lasers as a possible source for the electromagnetic 
radiation El. 



IV. QUANTUM SIGNATURES OF THE 
STABILIZED ASYMMETRIC STRETCH ORBIT 

After having shown the possibility to stabilize the 
asymmetric stretch orbit in the classical eZe configura- 
tion of helium, we discuss in this section to which ex- 
tent this classical stabilization manifests itself in the cor- 
responding quantum system. In principle, a full-blown 
three-dimensional treatment of doubly excited helium in 
the electromagnetic field would be desirable in this con- 
text. While such ab initio calculations can indeed be 
performed in the unperturbed two-electron atom (up to 
principal quantum numbers of the order of N ~ 15 ... 20 
|45| - which, as we shall see later on, should be sufficient 
for our purpose), they are still impractical in the pres- 
ence of a nonperturbative time-dependent electric field 
where the conservation of the total angular momentum 
is broken and the Floquet formalism is required to obtain 
wave packet states associated with the resonantly driven 
collinear orbit. Essential properties of these wave pack- 
ets, however, can be reproduced from a one-dimensional 
model of the two-electron atom which properly takes into 
account the dynamics along the field polarization axis. 
This is particularly the case if a magnetic field is ap- 
plied in order to provide stability in the transverse de- 
grees of freedom; the electrons are then strongly confined 
in the x—y plane and evolve according to an effective one- 
dimensional dynamics along the z axis. 

To correctly describe this dynamics, we demand that 
the one-dimensional model represents the exact quantum 
analog of the classical eZe configuration. This in partic- 
ular requires to take into account the full Coulomb inter- 
action between the charged particles. A smoothening of 
the Coulomb singularity, which is frequently employed in 
one-dimensional models of driven atoms (e.g., 46] ), may 
not be permitted here, since it would lead to a consid- 
erably different behaviour in the corresponding classical 
system. 

Consequently, we write the Hamiltonian that generates 
the quantum dynamics of the driven collinear configura- 
tion as 



H 



1 d 2 

2 9Ci 2 



2dd ~ Ci 



z 

~2 



1 



C1 + C2 



F . ( 8 d 
-— smut 5^-57- 
iuj \d(i dC,2 



(8) 



with Q\ = z\ and £2 = — zi the absolute values of the 
(Cartesian) coordinates of the electrons along the field 



polarization axis. The external field, parametrized by 
the amplitude F and the frequency w, is incorporated in 
the velocity gauge in order to ensure good convergence of 
the numerical calculation [47^. Effectively, the electrons 
appear here as distinguishable particles (with £ij C2 > 0), 
which is consistent with the classical impenetrability of 
the Coulomb singularity at the origin. 

Due to the temporal periodicity of the Hamiltonian, 
the Schrodinger problem represented by J3J is conve- 
niently treated in the framework of Floquet theory. 
The latter states that any solution V'tCCi^) of the 
Schrodinger equation can be written as 



d£C £ 4 £) (Ci,(2)e- 



-iS t 



(9) 



with time-independent complex expansion coefficients 
Cs, where £ are the quasienergies and tpf (C,x,Cs) — 
V , t+2 7 i7 ta ,(Ci> C2) t ne associated time-periodic quasienergy- 
eigenfunctions. The latter are determined by the time- 
independent eigenvalue equations 



(H + ku; - £)4 £) (Ci,C2) 



+ ^(C ) 1 (Ci,c 2 ) - C-i(Ci,c 2 )) 



(10) 



with 



H n = -- 



1 d 2 



1 d 2 



1 



2 ocl 2 acl 



d 



z z 

Ti ~ ~2 + C1TC2' 



2uj VdG d( 2 



which result from the Fourier series expansion 



# ) (Ci,C 2 ) 



k=— cx 



i kujt 



(ii) 



(12) 



For atomic systems, the ^-periodic Floquet spectrum 
of quasienergies £ is absolutely continuous: each bound 
state of the unperturbed atom is coupled to the atomic 
continuum via multiphoton transitions and therefore ap- 
pears as a resonance structure in the spectrum. In order 
to separate these spectral resonances from the flat back- 
ground, we employ the method of complex scaling (see 
48] for a recent review), which has proven its usefulness 
in a variety of numerical studies in atomic and molec- 
ular physics. Essentially, the electronic coordinates are 



complexified according to Q 1— » Qe l6 ', 



according to d/dQ 



~ W d/dQ (j 



and the momenta 
= 1,2). Formally, 



this scaling can be expressed by the (non-unitary) trans- 
formation 



# } (Ci,C 2 ) 



Wt (Cie ,C2e ) 



(13) 



of the quasienergy-eigenfunctions, as well as by the cor- 
responding transformation 



d d 
//(Cl ' C2 '9Ci"'9^' t 



(14) 
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of the Hamiltonian. As a result, we obtain a complex 
symmetric rather than Hermitian eigenvalue problem. 
Resonance structures in the continuous spectrum of the 
"real" , unrotatcd eigenvalue problem (6 = 0) appear now 
as discrete complex eigenvalues £ = E — iT/2, with the 
real and imaginary parts corresponding to the energies 
E and spectral widths T/2 (HWHM - the half width at 
half maximum) of the resonances, respectively. 

In perfect analogy to the quantum description of the 
driven Zee configuration of helium (see |2l| for more de- 
tails), the complex scaled Floquet Hamiltonian is now 
expanded in the product basis 

{^ ) (Ci)4 a) (C2):n,m>l} (15) 

composed of the real- valued Sturmian functions 

^(0- <±£f expUW (|), ,16) 
Vii a \ a J \ a J 

where the Ln denote the associated Laguerre polyno- 
mials. By the choice of this basis, the Hilbert space is 
effectively restricted to functions that scale at least lin- 
early with Q for small Q (j = 1> 2) and therefore do not 
exhibit divergent potential matrix elements. The scal- 
ing parameter a, which can be freely chosen in principle, 
determines the maximum spatial extension (as well as 
the minimum coarse graining) that is represented by a 
truncated basis set {s[ a \ . . . , S'jv' 1 }. Variations of a (as 
well as of the complex scaling angle 9) provide an effi- 
cient means to verify the numerical convergence of the 
calculation. 

In the unperturbed atomic system (F = 0), the eigen- 
states of H are characterized by a well-defined "parity" : 
they are either "even" or "odd" with respect to the ex- 
change of £i and (2- 

lKCi,C2) = ±^(Ca,Ci). (17) 

Within each of these two symmetry classes, the spectrum 
can be grouped into series that converge towards single 
ionization thresholds, which are labelled by the principal 
quantum number N of the inner electron. Naturally, this 
classification cannot be carried out in a perfectly rigorous 
way, since N is not a good quantum number. In practice, 
ambiguities arise as soon as states belonging to different 
series start to resonantly interact with each other, which 
is the case for N > 5 (see also [49^1. 

For each parity, the energetically lowest member of 
each series exhibits a pronounced localization in the 
vicinity of the classical asymmetric stretch orbit. This 
is illustrated in Fig. which shows the densities of the 
lowest even and odd states at principal quantum number 
N = 10. The probability distributions of these "asym- 
metric stretch states" are in fact rather similar to the pro- 
jections of three-dimensional intra-shell eigenfunctions 
(with maximum interelectronic angle) onto the £1 - £2 
subspace [25| , which are also localized along the collinear 
asymmetric stretch orbit. 
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FIG. 9: Probability densities of asymmetric stretch states in 
configuration space. Plotted are in (a) the energetically lowest 
"even" and in (b) the energetically lowest "odd" state at prin- 
cipal quantum number N = 10. We clearly see a localization 
along the classical asymmetric stretch orbit, which is shown 
by the underlying curve. With respect to Fig.^a), the orbit 
was scaled according to (j \— > (27rn llodc /S'o) 2 Cj (j> 2), with the 
number of nodes along the orbit given by n no( j c = 20 for the 
even and n no d c = 21 for the odd state. 



Due to the enforced node at £1 = £2, the energy E^- 
of the odd asymmetric stretch state is, at given N, con- 
siderably enhanced with respect to the energy E^+ of the 
even state, and lies approximately in the middle between 
En+ and £^+1)+- More precisely, the scaling laws (J2J 
predict the relation E ~ Eq{S/ Sq)~ 2 between the energy 
E and the action S (as defined by (J3J) in the semiclas- 
sical limit, where Eq ~ —1.50 and So ~ 18.8 are the 
respective values of the energy and action at frequency 
ijj = 1 (for the unperturbed asymmetric stretch orbit). 
Semiclassical quantization of the action along the orbit 
yields S ~ 2im no d c for the asymmetric stretch states. 
Here, n no d c denotes the number of nodes of the wave- 
function along the orbit (including the nodes at £1 = 
and C2 = 0), which is related to the principal quantum 
number N according to 

_ J 2N for even states , ■. 

n„odc - I 2 JV + 1 for odd states ' [ * } 

The resulting scaling 

£W±[a.u.] ~ £ (2™ nodc /S r 2 =i -13.4/n' ode (19) 

is, as shown in Fig. 1101 indeed encountered in the numer- 
ically calculated spectrum. 

For a resonant stabilization of the asymmetric stretch 
orbit associated with the state with node number n no de, 
the driving frequency needs to be set to u> = u~ where 
Vn naAx denotes the scaling parameter in J2J) that generates 
the transformation (J2J from the "scaled" orbit (at u> = 1) 
to the actual orbit underlying this state. Semiclassically, 
the action of the lowest quantized state within the stabi- 
lized asymmetric stretch island is given by 

S~27r(n nod0 + ( 7l + 72 )/2) (20) 
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FIG. 10: Energies of the asymmetric stretch states, as a func- 
tion of l/n^ ode with n no de the number of nodes along the 
asymmetric stretch orbit. The signs + and x denote "even" 
and "odd" states, respectively. The grey line shows the scal- 
ing E = — 13.4/n„ ode , which is expected in the semiclassical 
limit. 
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where 71 and 72 denote the winding numbers along the 
stabilized orbit (compare, e.g., with pin|). By virtue of 
5 = So^riaode) therefore, we infer 

^"„odc = 27r(n nodc + (71 + 7 2 )/2)/5 

~ 0.33(n nodo + 0.215) , (21) 

as evaluated for the numerically calculated values Sq ~ 
19.0, 71 ~ 0.174, and 72 — 0.256 (as obtained from the di- 
agonalization of the monodromy matrix) at Fo — 0.074. 
Comparison of the resulting approximate scaling uj ~ 
28/n^ odo with lfH?)l shows that w roughly corresponds to 
the level spacing between two consecutive asymmetric 
stretch states. This implies that these states are nearly 
degenerate in the resulting Floquet spectrum, and be- 
come strongly coupled to each other as soon as the time- 
periodic perturbation is switched on. Provided the size 
of the classical island is large enough, this coupling gives 
rise to a Floquet state that corresponds to a nondisper- 
sive two-electron wave packet anchored on the stabilized 
asymmetric stretch orbit. 

Contrary to the estimation based on the EBK crite- 
rion (see Section HIjl. such nondispersive wave packets 
are already found at principal quantum numbers N > 8. 
An example of such a wave packet, centered around the 
even asymmetric stretch state at N = 10, is shown in 
Fig. 1111 where the probability density of the correspond- 
ing Floquet state is plotted for the driving phases uit = 0, 
it/A, 7r/2, and 37r/4. The driving frequency was set to 
uj = 0.0034 a.u. ~ v^q, an( i the field amplitude was ac- 
cordingly chosen as F = 3.7 x 10~ 5 a.u. ~ V20F0, with 
Fo = 0.074 the scaled amplitude at which the size of the 
classical island is maximized. We clearly see that the 
wave packet faithfully traces the motion of the classical 
trajectory along the asymmetric stretch orbit. Significant 
contributions to the probability density far from the or- 
bit indicate, however, that the Floquet state is not yet 
perfectly localized on the island at this particular value 
of N. 



FIG. 11: Nondispersive two-electron wave packet anchored 
on the stabilized asymmetric stretch orbit. Plotted are the 
probability densities \ipt(Ci, C2 ) | 2 of the Floquet state ip t at 
(a) cut = 0, (b) at = tt/4, (c) cut = tt/2, (d) ut = 3tt/4. 
The driving frequency is set to uj — 0.0034 a.u. ~ v^o ( see 
Eq. ||2TJ), i.e., the wave packet is centered around the even 
asymmetric stretch state at N = 10. Accordingly, the field 
amplitude was chosen as F — 3.7 x 10 -5 a.u. ~ v^qFq with 
Jo = 0.074 (the underlying orbit was accordingly scaled, see 
Fig.|5Ja)). Significant contributions to the probability density 
far from the orbit reveal that the classical resonance island is, 
at this value of N, not yet big enough to support a fully 
localized quantum state. 

The spectral evolution from the unperturbed asym- 
metric stretch state to the nondispersive wave packet 
on the stabilized orbit is illustrated in Fig. I12f a) which 
shows the Floquet spectrum function of the field 
amplitude F. Clearly, the Floquet level associated with 
the stabilized asymmetric stretch orbit is shifted towards 
higher quasienergies with increasing F (see |2lj for the 
analogous case in Zee helium), which is in accordance 
with the fact that the energy of the classical orbit be- 
comes slightly enhanced as well. Note the occurence of 
huge avoided crossings (e.g., between the o and the □ 
states), which make it difficult to determine the states 
that contain significant contributions to the driven asym- 
metric stretch state (in practice, the states were identi- 
fied by a visual inspection of the probability density in 
configuration space). 

The corresponding ionization rates T/2 (HWHM) of 
these Floquet states are shown in Fig. I12f b). Despite 
appreciable admixtures of other, less stable components 
around F ~ 3.3 x 10~ 5 a.u. which "spoil" the ioniza- 
tion rate, it becomes apparent that a local minimum of 
T/2 occurs in the range 3 x 10~ 5 a.u. < F < 4 x 10~ 5 
a.u. where the cross section area of the classical island 
is maximized. At F — 3.7 x 10~ 5 a.u. in particular, the 
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FIG. 12: Evolution of the nondispersive wave packet state in 
the Floquet spectrum. The upper graph shows the quasiener- 
gies of the driven collinear atom at u = 0.0034 a.u. as a 
function of the field amplitude F. Plotted are the levels of 
resonances with ionization rate T/2 < 10 -4 that are coupled 
to the asymmetric stretch state by at most 10 photons. The 
marks o, □ and o indicate the Floquet states that share the 
contributions of the wave packet state associated with the sta- 
bilized asymmetric stretch orbit. The corresponding ioniza- 
tion rates (HWHM) of these states are shown in (b). Indeed, 
a local minimum of T/2, "spoiled" by near-resonant admix- 
tures of other, more unstable components, is encountered in 
the range F ~ 3 x 10 -5 . . . 4 x 10~ 5 a.u. where the size of the 
classical resonance island becomes maximal. 



ionization rate T/2 ~ 2.5 x 10~ 6 a.u. is obtained for the 
stabilized asymmetric stretch state. This implies that 
a wave packet that is initially prepared on the Floquet 
state shown in Fig. 1111 propagates, without spreading, 
along the asymmetric stretch orbit for about 5000 field 
cycles, until the population is completely distributed in 
the ionization continuum. 

As is clearly shown in Fig. I12f , the minimal ioniza- 
tion rate in the range 3 x 10~ 5 a.u. < F < 4 x 10~ 5 
a.u. is of the same order as the width of the unperturbed 
asymmetric stretch state. Globally, therefore, the exter- 
nal driving does not lead to a significant stabilization at 
N = 10, which can be traced back to the fact that the 
wave packet cannot be fully contained within the classical 
resonance island (see Fig. 1111) . A more pronounced sta- 
bilization effect is expected in the regime N ~ 20 . . . 30 
for which the EBK quantization criterion predicts the 
existence of fully localized quantum states on the res- 
onance island. Such quasi-bound states decay only via 
tunneling-induced couplings to the surrounding chaotic 
sea and to the ionization continuum, and should there- 
fore exhibit significantly lower ionization rates than the 
unperturbed asymmetric stretch states. Calculations in 
this regime, which require a substantially higher numer- 



ical effort than for N = 10, are presently under way. 



V. CONCLUSION 

In summary, we have shown that the collinear asym- 
metric stretch orbit of the classical helium atom can be 
stabilized by means of a resonant electromagnetic per- 
turbation. Within a reasonable range of field amplitudes 
(0.04 <Fo< 0.15 at u> = 1), the external driving induces 
a local regular island in the phase space of collinear mo- 
tion, in which the electrons perform stable quasiperiodic 
oscillations around the stabilized asymmetric stretch or- 
bit. The driven collinear atom is unstable with respect to 
deviations from collinearity. However, a static magnetic 
field, oriented parallel to the driving field polarization, 
can be used to stabilize the configuration by transver- 
sally confining the electrons to the vicinity of the field 
polarization axis. 

The quantum signatures of this classical stabiliza- 
tion phenomenon were elucidated by means of a one- 
dimensional model of the two-electron atom which cap- 
tures the essential properties of the dynamics along the 
field polarization axis. Within the Floquet spectrum of 
the driven atom, wave packet states were identified which 
follow, without spreading, the time-periodic Kepler-like 
motion along the stabilized asymmetric stretch orbit. 
Such nondispersive wave packets, which constitute the 
quantum analog of the classical nonlinear resonance, arise 
for symmetric double excitations with principal quantum 
number N > 8. At N = 10, a local maximum of the life 
time of the order of 5000 field cycles was found near the 
scaled field amplitude F ~ 0.07 where the cross section 
area of the classical resonance island is maximized. 

The experimental realization of such nondispersive 
wave packets is certainly difficult, but does not seem im- 
possible. Intrashcll "asymmetric stretch" -like states with 
maximum interelectronic angle were already populated 
in high-resolution photoionization experiments up to the 
principal quantum number N = 10 [3?! l3Sj . At that 
double excitation, the resonant driving of the asymmet- 
ric stretch orbit would correspond to a laser frequency of 
the order of u>/(2ir) ~ 10 13 Hz, while a laser intensity of 
the order of I ~ 10 8 W/cm 2 would be required for sta- 
bilization. On the other hand, the magnetic field that is 
necessary to enforce the transverse confinement along the 
field polarization axis would exceed 100 Tesla at N = 10, 
which is too large to be realized in table-top like exper- 
imental setups. It is not excluded, however, that wave 
packet states can be localized on the driven collinear or- 
bit even without the presence of an additional perturba- 
tion that stabilizes the dynamics in the transverse de- 
grees of freedom (in the same way as "scars" , which are 
anchored on unstable periodic orbits in chaotic systems 
H2). To examine to which extent such "scarred" wave 
packet states exist in the absence of a stabilizing mag- 
netic field (and to which extent the classically allowed 
decay via the transverse degrees of freedom would influ- 
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ence their life time), a more quantitative study, taking 
into account the dynamics within the two-dimensional 
configuration space, would be required. 
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APPENDIX A: NUMERICAL INTEGRATION OF 
THE CLASSICAL EQUATIONS OF MOTION 

In this Appendix, we describe the canonical 
"Kustaanheimo-Stiefel" transformation which is em- 
ployed in order to ensure stable numerical integration 
of Hamilton's equations of motion including electron- 
nucleus collisions. We restrict the description to the spe- 
cial case of collinear motion of the electrons in the eZe ar- 
rangement. Generalizations to two- or three-dimensional 
dynamics can be found in Refs. [S^ and [l^, respec- 
tively. 

In the subspace of collinear eZe motion, the Hamilto- 
nian of driven helium is given by 

Z Z 1 

+ \z\ + z 2 )F coscot 



H = 



Pli 



Zl 



Z-2 



zi - z 2 



. (A1) 

with zi > and z 2 < 0. Obviously, a direct numerical in- 
tegration of the equations of motion resulting from (|A1() 

leads to a diverging momentum p Z i ~ z- [i = 1,2) 
whenever one of the electrons hits the nucleus. We there- 
fore perform the canonical transformation 

Z\ h-> Qi = yfzj, Pzl | -> Pi = 2 y^Tp g i, 

Z 2 l-> Q 2 = V 22, Pz2 l-> Pi = 1y/-Z2Pz2, 

which is associated with the generating function 

FiPzl,Pz2,Ql,Q2) =PzlQl+Pz2Ql (A3) 

By construction, these new momentum variables Pi, P 2 
remain finite at electron-nucleus collisions. 

In addition to this canonical transformation, we intro- 
duce the new effective Hamiltonian 



(A2) 



H = H(Pi 1 P 2 ,Qi,Q2,E 1 t) 
= Q 2 1 Ql(H(p z i,p z2 ,z 1 ,z 2 ,t) 



-E) 

Qi) 



ziQiPt + QiPt) - z(Qi + Q z 2 ) - EQiQi 

QlQliQl -Qt)F coscut (A4) 



+ 



QiQ 2 



which defines, via the Hamiltonian equations of motion 
dPi dU 
dr dQ. k 
dQj = &H 
dr ^ dP l 



(A5) 
(A6) 



the time evolution of P%, P 2 , Q±, Q 2 with respect to the 
pscudo time r. Note that t and E appear as additional 
canonically conjugate variables, satisfying 



dE 
dt 



dt ~ QlQ2 dt 



m 

~8E 



(A7) 
(A8) 



Indeed, one can straightforwardly show that the com- 
bined set of differential equations l|A5HA8Jl leads to the 
same time evolution of Zi and pi as the original equations 
of motion derived from IjAlf) . provided the initial value 
of the energy variable is set to E — H\ t =o- Eqs. (|A5HA8p 
permit stable numerical integration over simple electron- 
nucleus collisions, and become singular only in the ex- 
ceptional event of a triple collision. 



APPENDIX B: IDENTIFICATION OF THE 
PERIODICALLY DRIVEN 
ASYMMETRIC-STRETCH ORBIT 

The periodic orbits of the driven eZe configuration of 
helium are identified using the so-called stability trans- 
formation (ST) method [£l EB, E3 ■ The latter al- 
lows to detect unstable periodic orbits of length p in the 
n-dimcnsional chaotic time-discrete systems 



f : X i — > X' = f (X) 



(Bl) 



i.e. the fixed points of the p times iterated map X i — > 
fW(X). 

To this end, a stability transformed system s is intro- 
duced according to 



s : X i — > s(X) = X + AC[f (p) (X) - X] , 



(B2) 



where < A <C 1 is a scalar parameter and C is a con- 
stant, regular, and real matrix. While the positions of 
the fixed points of s are exactly the same as in the origi- 
nal system f( p \ their stability properties depend on the 
matrix C and the parameter A. The goal is to choose 
these parameters in such a way that the fixed point Xo 
to be located becomes stable with respect to the map 
s; the more unstable Xo is in the original system, the 
smaller A has to be chosen in order to achieve stabiliza- 
tion jEJ 0, HE X is then obtained as a point 
of convergence when propagating a properly chosen set 
of initial points with the stability transformed system s. 
The most dominating advantages of the algorithm are the 
large extensions of the basins of attraction of the individ- 
ual fixed points and the fast convergence of trajectories 
far from these fixed points. 

In general, a single transformation s is not sufficient to 
stabilize all fixed points of the original system. The above 
procedure is therefore applied to a set of transforma- 
tions Si with different matrices C, (i = 1,2, N). For 
our particular case of motion within the four-dimensional 
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phase space spanned by (zi, Z2,p z \,Pz2), we find that the 
set of matrices 




-1 



(B3) 

is sufficient for 



with "•" = 0, "+" = +1, and "-" 
the identification of all fixed points. 

For time-continuous systems, periodic orbits are repre- 
sented by fixed points of a suitably defined Poincare map 
f (X), which permits their detection by the ST method 
|57|. In our case of a periodically driven system with 
two degrees of freedom, the Poincare map is most conve- 
niently chosen as the stroboscopic map that transforms 
the phase space variables to their propagated values after 
one period T — 2tt/lo of the driving: 



f : (zi,z 2 ,p z i,Pz2)\t ' — ► (zi,z 2 ,p z -L,p z2 )\t+i 



(B4) 



An initial distribution that is chosen uniform in phase 
space (with the additional limitation — 5 < E < for the 



initial energy E) proved to be favourable for the detection 
of large sets of periodic orbits with unlimited lengths. 

For polishing-up the coordinates X = {z\, Z2-,p z i,p Z 2) 
of the periodic orbits, a Newton algorithm is applied af- 
ter convergence of the ST method up to a certain ac- 
curacy e ~ 10~ 2 . To this end, the monodromy ma- 
trix associated with the actual fixed point Xo is ap- 
proximately calculated by propagating the four points 
Xi = X + rjei (i = 1, ...,4) under the Poincare map f, 
where are orthogonal unit vectors in phase space and 
< 77 <C 1. The individual columns of the mon- 
odromy matrix M = (mi, ni2, 1x13, 1114) are then given 
by rrii = (f(Xi) — f(X))/?7. This monodromy matrix is 
then used to improve the accuracy of the fixed point by 
the Newton method. 

In a similar way, the Newton method can also be em- 
ployed to identify a given periodic orbit at the field am- 
plitude F + SF, once the corresponding orbit at field 
amplitude F is converged. In that case, the fixed point 
Xq at F would be used as starting point for the search 
of the fixed point 'X^f +SF ' at F + 5F (a robust damping 
of the Newton algorithm by a factor p « 0.1 is neces- 
sary there to ensure stable convergence H^)- This proce- 
dure turned out to be particularly efficient for identifying 
(within a single run) the asymmetric stretch orbits at all 
field amplitudes within the range < F < 1 (a stepsize 
SF = 0.001 was used in practice). Also other periodic or- 
bits, which were found by the ST method, were "traced" 
in this way for increasing field amplitude F. While their 
Lyapunov exponents may exhibit local minima at finite 
F, none of them - except the asymmetric stretch orbit - 
was found to become stable. 
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